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Abstract Electronic excitation energy transfer along 
a molecular wire depends on the relative orientation of 
the electronic transition dipole moments of neighbor- 
ing chromophores. In porphycenes this orientation is 
changed upon double proton transfer in the electronic 
ground state. We explore the possibility to trigger such 
a double proton transfer reaction by means of an in- 
frared pump-dump laser control scheme. To this end 
a quantum chemical characterization of an asymmetri- 
cally substituted porphycene is performed using density 
functional theory. Ground state geometries, the topol- 
ogy of the potential energy surface for double proton 
transfer, and So— >Si transition energies are compared 
with the parent compound porphycene and a symmetric 
derivative. Employing a simple two-dimensional model 
for the double proton transfer, which incorporates se- 
quential and concerted motions, quantum dynamics sim- 
ulations of the laser driven dynamics are performed 
which demonstrate tautomerization control. Based on 
the orientation of the transition dipole moments this 
tautomerization may lead to an estimated change in 
the Forster transfer coupling of about 60%. 
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1 Introduction 

Molecular photonic wires are of considerable interest 
as biomimetic models for photosynthetic light harvest- 
ing systems and as molecular devices [T]. Lindsey et. 
al. have been the first to demonstrate directed energy 
transfer through a porphyrin-based wire after selective 
excitation of a boron-dipyrrin donor 2 . Subsequently, 
they described a number of different systems on the 
basis of multiporphyrin arrays [HHUE] including, for in- 
stance, shape persistent cyclic architectures [6] or op- 
toelectronic gates 7,8 (for a review, see also Ref. [9]). 
The energy transfer in these structures is dominated by 
a through-bond mechanism with some dependence on 
the nature of the bridge linking the chromophores |10) . 
A different design strategy has been followed by Sauer 
and coworkers who have used the backbone of double- 
stranded DNA as a scaffold for positioning highly fluo- 
rescent chromophores [TTl[T2l[T5irLl] (see also Ref. |T5] ) . 
Here, Forster type dipole-dipole interaction driven trans- 
fer is the responsible mechanism for energy transfer 
with efficiencies close to 100% over a range of ~14 nm. 

The utilization of molecular wires in optoelectronic 
devices requires to have at hand a means for switching 
the energy transfer dynamics. Various methods such as 
optical excitation using UV radiation |16) , electrochem- 
ical oxidation [8l lT71[T8l or changing the pH value of the 
medium [T5] have been suggested (for an overview, see 
also Refs. [Tll20|h This includes even the realization of 
photochemical logic gates as reported in Ref. Re- 
cently, we have proposed another route to switch exci- 
tation energy transfer which could operate in the ultra- 
fast regime of a few picoseconds [22]. The idea builds 
on the observation that the direction of the electronic 
transition dipole moment, which is responsible for the 
coupling between chromophores in the Forster mecha- 
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nism [23] , depends on the positions of the two H atoms 
in the porphyrin like building blocks of molecular wires. 
Thus switching between the ground state tautomers, i.e. 
triggering the double proton transfer (DPT) by means 
of an appropriately designed laser field in the infrared 
(IR) domain enables one to modify the strength of the 
Forster transfer coupling. 

The aim of the present contribution is to extend 
the model study of Ref. |H] towards a specific molec- 
ular system and to show exemplarily to what extent 
laser-triggered DPT can influence the Forster transfer 
coupling. Specifically, we will focus on porphycene like 
structures, which have rather strong hydrogen bonds 
(HBs) whose properties are widely adjustable by the 
structural design [24]. For the laser-driven DPT the 
shape of the potential energy landscape is rather im- 
portant. Besides energy barriers which need to be sur- 
mounted by excitation with a few IR photons, reactants 
and products should be energetically different to pre- 
vent an efficient backreaction via tunneling. Asymmet- 
rically substituted porphycenes seem to meet these re- 
quirements. Below we will consider 9-actetoxy-2,7,12,I7- 
tetra-n-propylporphycene (3) whose spectroscopy had 
been studied by Waluk and coworkers [25]- Thus, be- 
sides the specific goal of demonstrating a laser-driven 
switch for molecular photonic wires, our quantum chem- 
ical studies of the ground state properties and the Sq— >S± 
excitation spectrum will shed some light on the effect 
of symmetric (1) and asymmetric (3) substitution of 
the porphycene parent compound (1) and hence will 
provide the theoretical basis for the assignment of the 
absorption spectrum. 

The paper is organized as follows: Section [5] gives a 
brief summary of the quantum chemical and quantum 
dynamical methods. The results Section [3] starts with 
the discussion of the geometries and energetics of the 
stationary points of porphycene and its derivatives (see 
Fig. [I}. Based on this information a two-dimensional 
(2D) model Hamiltonian is constructed in Section 13.21 
which comprises coordinates of concerted and sequen- 
tial DPT. Subsequently, the So— ^Si excitation spectrum 
is discussed in terms of the involved molecular orbitals 
and transition dipole moments. The laser driven DPT 
dynamics will be presented in Section 13.51 and the final 
Section [4] concludes with a summary. 



2 Methods 




Fig. 1 Reference structures of considered porphycenes: 1 free 
base porphycene, 2 tetra-n-propyleporphycene, and 3 more stable 
trans form (Trl) of tetra-n-propyleacetoxyporphycenc. 

validated by means of frequency calculations. Excita- 
tion energies are obtained from time-dependent DFT 
(TDDFT) calculations at the optimized ground state 
geometries (vertical excitation) and using the same func- 
tional and basis set. All quantum chemical calculations 
are performed with the Gaussian 03 suite of programs 

The two-dimensional time-dependent Schrodinger e- 
quation has been solved using the multiconfiguration 
time-dependent Hartree (MCTDH) approach 27, 28] as 
implemented in the Heidelberg program package 29 . 
The grid was chosen according to the harmonic oscil- 
lator discrete variable representation (64 points within 
[-2.5:2.5]ao). The actual propagation is performed using 
the variable mean field scheme in combination with a 
6-th order Adams-Bashforth-Moulton integrator. Using 
20 single particle functions per degree of freedom the 
largest natural orbital populations have been typically 
on the order of I0~ 6 . Selected cigenstatcs of the time- 
independent Hamiltonian, ifi , are obtained by improved 
relaxation [30 . 

3 Results 



Geometry optimization for the electronic ground state 3.1 Optimized geometries 
has been performed using density functional theory (DFT) 

employing the B3LYP hybrid functional together with Ground state optimized structures of 2 show that this 

a 6-31+G(d,p) basis set. Stationary points have been molecule like 1 is planar in accord with Ref. [23]. Asym- 
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Fig. 2 Optimized geometries (B3LYP/6-31+G(d,p)) of all sta- 
tionary points of 3 which are relevant for DPT. 

metric substitution, however, causes 3 to become non- 
planar with respect to the acetoxy group. Overall 3 
has two global trans minima, Trl and Tr2, two local 
cis minima, Cisl and Cis2, four local maxima which 
correspond to first order saddle points, Tsll to Ts22, 
and one global maximum being a second order saddle 
point, SP2. All these structures are shown in Fig.[2]and 
the energetics is compiled in Tab. [TJ In terms of the 
double proton transfer, SP2 is passed in the concerted 
mechanism, while pathways like Trl-Tsll-Cisl-Ts2I- 
Tr2 are followed during sequential transfer. Finally, we 
note from the dipole moments given in Tab. Q] that 3 is 
considerably more polar than 2. 

The energetic difference between trans and cis, which 
is 804 cm -1 for 1, is largely unaffected by the symmet- 
ric substitution in 2. The same holds true for the sec- 
ond order transition state SP2, i.e. the barrier height 
for concerted DPT is about 2260 cm -1 . The barrier for 
stepwise proton transfer which is 1637 cm _1 for 1 de- 
creases slightly to 1514 cm _1 in 2. Upon asymmetric 
substitution with an acetoxy group in the 9 position in 
3 most notably in comparison with 2 is the asymmetry 
of the two trans tautomers which amounts to 38 cm -1 . 
This rather small value is in good accord with the ex- 
perimental finding of indistinguishable tautomers [25] . 
Further, the SP2 transition state is lowered by about 
4% and the first order saddle points and cis minima are 
changed in a way such as to make a stepwise DPT via 
Cisl energetically more preferable. 



Table 1 Energies and absolute value of dipole moments of sta- 
tionary points for the tautomerism of the molecules shown in Fig. 
□ at the B3LYP/6-31+G(d,p) level of theory (values for 1 from 
Ref. 131)) ■ 



molecule 


structure 


AE (cm 1 ) 


dipole (Debye) 


3 


Trl 





1.633 




Tr2 


38 


1.481 




Cisl 


726 


2.042 




Cis2 


838 


2.218 




SP2 


2169 


1.635 




Tsll 


1329 


1.539 




Tsl2 


1505 


1.865 




Ts21 


1479 


1.690 




Ts22 


1461 


1.834 


2 


Tr 





0.002 




Cis 


817 


1.452 




SP2 


2263 


0.000 




SP1 


1514 


0.770 


1 


Tr 









Cis 


804 






SP2 


2256 






SP1 


1637 





The results on the HB geometries of 3 are compiled 
in Tab. [2] which also contain values for the trans tau- 
tomers of 1 and 2. Comparing the most stable trans 
forms of the three molecules the differences in HB pa- 
rameters are rather small. Most notably again is the 
asymmetry, i.e. in Trl the HB next to the acetoxy sub- 
stitution site (N1-H3-N2) is slightly shorter and sub- 
stantially more nonlinear as compared with the N3-H!,- 
N4 HB. The geometric signature of the energetically 
more preferable path for stepwise DPT can be found in 
the stronger contraction of the N3-Hf,-N4 HB in Tsll as 
compared with the Ni-H Q -N 2 HB in Tsl2. The transi- 
tion states for the second transfer step, Ts21 and Ts22, 
on the other hand, have almost the same HB parame- 
ters. 



3.2 Ground state model Hamiltonian 

The potential energy surface for DPT in the electronic 
ground state will be modeled using a simple Hamilto- 
nian which was introduced by Smedarchina and cowork- 
ers [331134] and which has been supplemented by asym- 
metric terms to mimic systems like 3 in Ref. |22j . As- 
suming that only the proton coordinates are active on 
a surface where all other coordinates are adiabatically 
adjusted we have the following Hamiltonian (m = 2mjj) 

2m { dxl + dxl ) + Usym [Xs ' Xa) 

"f" ^asym(^-S5 -^a) • (1) 

Here, symmetric and asymmetric transfer coordinates 
x s and x a , respectively, have been introduced. They are 




1 



Table 2 Geometries of the HBs in 3 (and for Tr in 1 and 2) at the B3LYP/6-31+G(d,p) stationary points (bond lengths in Angstrom, 
angles in degrees, data for 1 are taken from Ref. [32 J. 



structure 


Ni-Ha 


N 2 -H a 


Ni-N 2 


ZNi-H a -N2 


N 3 -H 6 


N 4 -H 6 


N3-N4 


ZN 3 -H b -N 4 


Trl 


1.05 


1.68 


2.65 


152.9 


1.05 


1.70 


2.67 


150.9 


Tr2 


1.68 


1.05 


2.65 


153.0 


1.67 


1.05 


2.66 


154.9 


Cisl 


1.61 


1.07 


2.61 


155.1 


1.06 


1.64 


2.62 


152.6 


Cis2 


1.07 


1.61 


2.61 


154.0 


1.60 


1.07 


2.62 


156.5 


SP2 


1.30 


1.26 


2.51 


159.4 


1.28 


1.27 


2.50 


159.1 


Tsll 


1.06 


1.61 


2.60 


153.1 


1.26 


1.30 


2.52 


159.2 


Tsl2 


1.30 


1.26 


2.52 


159.6 


1.57 


1.07 


2.58 


154.8 


Ts21 


1.26 


1.31 


2.52 


159.2 


1.07 


1.57 


2.58 


155.3 


Ts22 


1.57 


1.07 


2.59 


156.8 


1.31 


1.25 


2.52 


158.9 


2 


1.05 


1.69 


2.67 


152.6 










1 


1.05 


1.68 


2.66 


152.8 











defined with respect the to linear displacements of the 
two H atoms along the hydrogen bonds, x\ and 22, as 
follows: x s — (x\ + X2)l2 and x a = (xi — X2)j2. The 
symmetric potential, U sym (x s , x a ), originates from two 
bilinearly coupled quartic oscillators when expressed in 
the coordinates of the individual H atoms X\ and X2 ■ 
After introduction of the transfer coordinates x s and 
x a this gives: 

U sym {x s ,x a ) = 2U + -4 [(9 - *)x 2 a ~ (9 + 4)^] 
x Q 

+ 2 -^« + < + 6*2*2). (2) 

x 

Here g is the coupling constant, ±Xq are the minima for 
the uncoupled double minimum potential which has a 
barrier height of Uq. Asymmetry can be introduced by 
the following term: 

Uo a c i s U 

t^asym v*s ? *a ) — *s T X a (oj 

Xq Xq 

where atrans and a c i s are dimensionless parameters char- 
acterizing the detuning between the trans and cis states, 
respectively. 

The parameters for this Hamiltonian have been ob- 
tained by fitting to the quantum chemical results of 
Tab. [U The resulting potential energy surface is plot- 
ted in Fig. [3] and the fitting parameters are given in 
the figure caption. Note that while the minima and the 
SP2 maximum arc in good agreement with the quan- 
tum chemistry results, the simple form of the potential 
does not provide enough flexibility for obtaining all first 
order transition states at the same time with compara- 
ble accuracy. The deviations range from 51 cm -1 to 169 
cm , which is, however, still acceptable given the over- 
all approximate nature of the treatment. 

For the interaction with the laser field we assume 
that the permanent dipole moment depends only lin- 
early on the coordinates, i.e., 

fffieldO) = -(d a x a + d s x s )E{t) (4) 



Table 3 Excitation energies (in cm 1 ) and oscillator strengths 
(in parenthesis) obtained at the TD-DFT B3LYP/6-31+G(d,p) 
level of theory. The calculated values for the different structures 
are compared with the experimental assignment from Refs. 1251 
(2,3) and |37| (1). Note that the experimental values correspond 
to the respective 0-0 transition. 



struct. 


calc. exp. 


leading excitation 


3(Trl) 
3(Tr2) 
2 
1 


17392 (0.15) 15708 
17670 (0.13) 15993 
17536 (0.14) 15983 
17904 (0.12) 16000 


144-5-145(0.58) 
143-5-145(0.58) 
128-5-130(0.48) 
80-5-82(0.56) 



where d a / s is the derivative of the dipole moment with 
respect to x a i s . The permanent dipole moment changes 
mostly along one polarization direction (x, cf. Fig. [SJ. 
Assuming the laser being linearly polarized along that 
direction and making a linear interpolation between the 
minimum configurations we obtain d s — — 0.028e and 
d a = 0.150e. 

For the laser field we will use the following form 
E(t) = E o,i cosM) exp(-(f - t 0ii ) 2 /2 T f) (5) 

i=l,2 

which facilitates the consideration of the so-called pump- 
dump control mechanism |35J[36J. In Eq. ([5]) Eq^ is the 
amplitude, a;, the carrier frequency, and the temporal 
width of the pulse centered at io,i- 

3.3 Molecular orbitals and excited states 

In the following we will discuss the So— >Si transition, 
which contributes to the Q absorption band. Results of 
TDDFT calculations are compared with experimental 
data in Tab. 3 for the different compounds. First, we no- 
tice that TDDFT gives excitation energies which exceed 
the experimental ones by about 1500 - 2100 cm -1 . Most 
of this difference must be attributed to the method it- 
self (see, e.g., systematic study in Ref. [38 ). Further 
deviations might be due to the fact that the experi- 
mental values correspond to the 0-0 transitions whereas 
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Fig. 4 Molecular orbitals relevant for the So— >Si transition (cf. Tab. 3) for molecules 1, 2, and the two tautomers of 3. The orbital 
numbers are 80-83 (1), 128-131 (2), and 143-146 (3). 



the calculations are performed for a vertical excitation. 
Comparing the experimental spectra of porphycences 
in gas and condensed phases the environmental effects 
(e.g., a nitrogen matrix in Ref. [25] ) is negligible [24] . 

Let us focus on the tautomerization of 3. As dis- 
cussed in Section 13.11 the energetic asymmetry in the 
electronic ground state is rather small. The measured 
splitting in the So— >Si transition, however, gave evi- 
dence for a considerable asymmetry in the excited states 



of 285 cm -1 [23J. The assignment of the most stable 
tautomer in the Si excited state has been based on a 
semiempirical ZINDO/S calculation. Considering Tab. 
3 we notice that the agreement of the splitting in the 
So— >S\ transitions (278 cm -1 ) is rather good thus giv- 
ing strong support for the experimental conclusions in 
Ref. [23]. 

The So— t-Si transitions are of tt — » ir* type with par- 
ticipation of the 7T (HOMO, HOMO-1) and tt* (LUMO, 



Fig. 3 Two-dimensional model potential according to Eq. Q for 
the DPT in 3. The parameters are Uo = 876 cm -1 , xq = 0.624ao, 
g = 0.435, atrans = 0.0095, and a cis = 0.034. The barrier height 
for concerted transfer is 2169 cm -1 . The energetic difference be- 
tween the two trans minima is 38 cm -1 and the energies of the cis 
minima are 725 cm - 1 and 838 cm" 1 (contour lines are in steps 
of 400 cm -1 from 400 cm -1 to 8800 cm -1 ). 



LUMO+1) molecular orbitals (MOs). These orbitals 
are shown for the different molecules in Fig. @] and 
orbital coefficients are compiled in Tab. 3. The main 
contributions to the So— >Si transition in 1 and 2 is of 
HOMO-f -> LUMO type (see also Ref. [39]). The same 
holds true for the less stable tautomer Tr2 of 3. In- 
terestingly, in tautomer Trl of 3 this transition is of 
HOMO ->■ LUMO type. Comparison of the MOs of Trl 
and Tr2 shows that the order of HOMO and HOMO-1 
is reversed upon tautomerization. 

A closer look at the orbitals reveals the effect of 
symmetric and asymmetric substitutions. Going from 
1 to 2 some electron density appears at those alkyl 
chains which are opposite to the position of the hy- 
drogen bonded H atom. Overall, however, the MOs are 
rather similar, i.e. in the HOMO there is a p-orbital at 
those nitrogens, which are bonded to the central hydro- 
gens whereas there is no electron density at the "free" 
nitrogen atoms. In the HOMO-I this situation is re- 
versed. Asymmetric substitution causes a disturbance 
of the porphycene 7r-system by the acetoxy group with 
the effect being different for the two tautomers. Tr2 of 
3 resembles the HOMO in 1 and 2 with some electron 
density at the acetoxy group. In Trl of 3 the presence 
of the latter electron density renders the Tr2-HOMO-I 
to become more stable than the Tr2-HOMO. 




Fig. 5 Transition dipole moment vectors for structures l(left) 
(1.252,-0.950), 2(left) (-1.489,0.672), 3(Trl) (-1.652,0.356), and 
3(Tr2) (-1.042, 1.174) (in Debye). 

3.4 Transition dipole orientation 

The directions of the transition dipole moments for the 
So— >Si transitions of the different molecules are com- 
piled in Fig. [5] Comparing 1 and 2 we notice that sym- 
metric substitution has almost no effect on the direction 
of the transition dipole moments. Furthermore, both 
cases show the well-known 90 degree change of the tran- 
sition dipole moment upon tautomerization. The pres- 
ence of the acetoxy group in 3, however, changes this 
picture due to the polarity introduced by the oxygen 
atoms. Most notably is the fact that the tautomeriza- 
tion causes a tilt of the transition dipole moment by 
about 50 degrees only. 

Let us assume that 3 is built into a molecular pho- 
tonic wire. The Forster energy transfer coupling be- 
tween different molecular subunits of the wire called 
donor (D) and acceptor (A) is given by [23] 

Jda = K DA d[° ) 4o } ■ (6) 

n DA 

Here, d^ D ^ is the magnitude of the transition dipole 
moment for the So— S"Si transition and Rda is the dis- 
tance between donor and acceptor. Most important for 
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the present discussion is the orientational factor given 

by 

k da = n A ■ n B - 3(e DA ■ n D )(e DA ■ n A ) , (7) 

where n D / A and ejj A are the unit vectors pointing in 
the directions of the transition dipole moment of D/A 
and the center to center DA distance, respectively. 

In order to estimate the effect of DPT on the Forster 
transfer coupling and in particular on the orientational 
factor we suppose that 3 is built into the wire such that 
the neighboring unit has a transition dipole moment 
like structure 1. This could be the case, for instance, 
for the widely used porphyrins. Let us further assume 
that ejj A points along the x-direction. In this case we 
obtain K DA (Trl)=1.43 and Kz 5j4 (Tr2)=0.60. In other 
words the transfer coupling changes by about 60% as a 
consequence of the tautomerization. 
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Fig. 6 Energy level scheme of the model Hamiltonian for 3 in- 
cluding the probability densities of those states which dominate 
the dynamics. The lengths of the horizontal lines indicate the 
area of derealization of the wave function along the x a coordi- 
nate. The green line corresponds to the energy of the SP2 bar- 
rier. The vibrational eigenstates have been obtained using the 
improved relaxation scheme within MCTDH I3UI . The pathway 
of laser-driven pump-dump control is marked by vertical arrows. 
The density plots cover the range [— 1 : l]ao along the vertical x a 
and the horizontal x s axes. 



3.5 Laser control of DPT in the electronic ground state 

In the following we will investigate the possibility to 
control the DPT in 3 using a simple two pulse pump- 
dump scheme involving an intermediate state which is 
energetically above the reaction barrier for concerted 
DPT (cf. Fig. [5]). The initial state is localized mostly 
in the Trl well, while the final state lies mostly in the 
Tr2 well. Due to the low barrier and the small asymme- 
try both states have finite probabilities in the respec- 
tive other wells. The goal is to find an electric field of 
the form given by Eq. J5]), for which the first (pump) 



pulse populates the intermediate state and the second 
(dump) pulse triggers a transition from this intermedi- 
ate to the reactant state. This scheme operates in com- 
plete analogy to other isomerization reactions studied 
previously, such as the single H-atom transfer in mal- 
onaldehyde |36] or the cope rearrangement in semibul- 
lvalenes [35] ■ The parameters for the present case of 
DPT were chosen as follows: First, the field strength is 
fixed at a certain value and the pulse length is increased 
starting from 800 fs until the intermediate state popula- 
tion reaches a maximum. For this pulse length the field 
amplitude is changed around its starting value to check 
whether the intermediate state population can be in- 
creased. The carrier frequency is chosen to correspond 
to the respective transition frequency. The same proce- 
dure is followed independently for the dump pulse by 
letting it start immediately after the pump pulse. This 
gave a total pulse duration of 5 ps during which the 
population of the intermediate state showed a plateau, 
which allowed us to shift the dump pulse such that 
pump and dump pulse partly overlap (cf. Fig. [7]). 

The progress of the laser-driven reaction can be fol- 
lowed by calculating the populations of vibrational eigen- 
states of the model Hamiltonian. A representation of 
the energy level scheme is given in Fig. [5] which also 
shows the probability density of those states which dom- 
inate the dynamics. The chosen above-barrier transition 
state involves an excitation of the asymmetric stretch- 
ing coordinate x a which triggers the reaction to oc- 
cur mostly along the stepwise pathway (see also Ref. 
[22]). Furthermore, it is seen that the probability den- 
sity has a higher weight in that part of the potential 
where x a < 0, implying that there is a preference for 
the Cisf minimum to be visited during the dynamics. 
According to our previous study [32] this is due to the 
higher stability of Cisf as compared with Cis2. Finally, 
we note that there are several states in the region of the 
chosen intermediate state (p4, thus state selective pop- 
ulation will be more easily achieved using long pulses 
as compared, e.g., with the NH- vibrational stretching 
frequency. 

The optimized field and selected state populations 
are shown in Fig. [7] The pump pulse excites the system 
from the ground state ground state, ipo, to the inter- 
mediate state tp4, while the dump pulse de-excites this 
population to the target state tpi . The maximum inter- 
mediate state population is 91 %, reached around 3.4 
ps. At this time, however, the dump pulse is already 
transferring population from 994 to the target state ipi. 
The overall duration of the laser-induced DPT process 
is 4.5 ps and the final population of the target state is 
95 %. Two minor features seen in Fig. [7J deserve a com- 
ment. Around 2.4 ps there is some transient population 
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of the target state already due to the pump pulse. Later 
on around 3.6 ps the dump pulse also is seen to further 
depopulate the initial ground state. Both effects can 
be traced to the fact that the energetic separation be- 
tween the two transitions is only 56 cm -1 , i.e. they can 
be simultaneously addressed by the laser pulse spectra. 
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Fig. 7 Upper panel: Electric field for the pump-dump control 
scheme (pulse parameters: Eo,i = 7 mEjj/eao), uii = 1111 cm -1 , 
t ,i = 2200 fs, and t\ = 650 fs; _E ,2 = 4 mE h /ea ), uj 2 = 1055 
cm -1 , to,2 = 3750 fs, and ri = 200 fs). Lower panel: Population 
dynamics of the states <po, ipi, and 1^4 



4 Conclusions 

In summary, we have investigated the possibility to use 
IR laser-controlled double proton transfer in the elec- 
tronic ground state of asymmetrically substituted por- 
phycene as a way of influencing the direction of the 
So— ^Si transition dipole moment. This gives a means 
for changing the dipole-dipole interaction between neigh- 
boring chromophores in a molecular photonic wire. The 
estimated effect amounts to a 60% change of the trans- 
fer coupling which could be used to design a molecular 
switch operating on a time scale of a few picoseconds 
and being, in principle, reversible. 

Needless to say that this proof-of-principle study 
calls for several extensions. First, the two-dimensional 
model of double proton transfer needs to be extended 



to account for heavy atom motions of the scaffold, e.g., 
by using the Cartesian reaction surface approach (40j . 
Indeed Waluk and coworkers showed that the coupling 
to this type of motions may have a strong influence on 
the dynamics [4T1I421I43] . It is to be expected that in- 
tramolecular vibrational energy redistribution is effec- 
tive on the time scale of a few picoseconds covered in 
the present study, i.e. either substantially shorter pulses 
have to be used or more sophisticated pulse forms need 
to be devised, which, e.g., make use of the promoting 
character of certain scaffold modes. Second, a suitable 
design for a molecular photonic wire containing a DPT 
switching unit has to be found. This implies not only a 
favorable orientation of the transition dipoles of neigh- 
boring units but also the dominance of through space 
Forster transfer over through bond energy transfer. 
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